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ABSTRACT 

We compute the polarization of the Lja line photons emerging from an 
anisotropically expanding and optically thick medium, which is expected 
to operate in many Lja emitting objects including the primeval galaxy 
DLA 2233+131 and Lyman break galaxies. In the case of a highly optically 
thick medium, the escape of resonance line photons is achieved by a large 
number of resonant local scatterings followed by a small number of scatterings 
in the damping wing. We show that some polarization can develop because the 
wing scatterings are coupled with strong spatial diffusion which depends on the 
scattering geometry and kinematics. The case of a slab with a finite scattering 
optical depth and expansion velocity of ~ 100 km s" 1 is investigated and it is 
found that Lja photons are emergent with the linear degree of polarization up 
to 10 per cent when the typical scattering optical depth r^lO 5 . We subsequently 
investigate the polarization of Lya photons emerging from a spherical shell 
obscured partially by an opaque component and we obtain ~ 5 per cent of 
polarization. It is proposed that a positive detection of polarized Lja with 
P-Cygni type profile from cosmological objects can be a strong test of the 
expanding shell structure obscured by a disk-like component. 

Subject headings: polarization : radiative transfer - cosmology : galaxy : 
individual DLA 2233+131 
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1. Introduction 

Various astronomical objects in the cosmological scales show P-Cygni type profiles in 
the Lya emission. These objects include the most remote galaxy at z — 4.92 gravitationally 
lensed by CL1358+62 (Franx et al. 1997), high z galaxies observed with the Hubble Space 
Telescope and the Keck telescopes (Steidel et al. 1996, Giavalisco et al. 1996, Lowenthal 
et al. 1997) and the damped Lyman a (hereafter DLA) candidates (Djorgovski et al. 
1996, 1997). Similar P-Cygni Lya profiles are found in nearby starburst galaxies, which are 
sometimes classified as Wolf-Rayet galaxies, blue compact galaxies, or H II galaxies (e.g. 
Kunth et al. 1996, Heckman and Leitherer 1997, Sahu and Blades 1997, Lequeux et al. 
1995, Legrand et al. 1997, etc.). The column density N HI of neutral hydrogen in these 
systems is usually found to be in the range N HI ~ io 19 ~ 21 cm~ 2 . 

The primeval galaxies or the first star clusters expected to be found at z > 5 epoch 
may possess a central super star cluster surrounded by neutral hydrogen of high column 
density (Haiman and Loeb 1997a,b). These surrounding layers can be accelerated by the 
expanding H II region just outside the super star cluster. It is hoped that in the near future 
with the advent of the Next Generation Space Telescope (NGST), the infrared spectra 
of these infant galaxies will be accessible and that the observational confirmation of the 
ubiquity of P-Cygni type Lja profiles may test the above hypothesis. 

Ahn and Lee (1998, hereafter AL98) investigated the Lya line formation in a thick 
and expanding medium. It was emphasized that the profile formation should be studied 
by accurately computing the contributions from photons back scattered by receding 
medium and wing-scattered photons (see also Legrand et al. 1997). It is well known 
that the properties of the Lya photons scattered in the damping wing are characterized 
by the Rayleigh phase function (e.g. Stenflo 1980). This is in contrast with the degree 
of polarization p = resulting from a resonance transition between 1Si/ 2 and 2P 1 / 2 and 
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p = 3/7 obtained for the 15*1/2 and 2P 3 / 2 transition (Lee et al. 1994). 

In a moderately thick and static medium a negligibly polarized flux is expected because 
the photons are locally scattered many times and get isotropized before they escape to the 
observer (e.g. Lee 1994). However, in a very thick medium, the escape is achieved by a 
large number of local resonant scatterings followed by a small number of scatterings in the 
damping wing. Hence, in the wing regime the spatial diffusion becomes important and the 
radiation field may get anisotropic depending on the scattering geometry. Therefore, the 
emergent line photons may get polarized and also anisotropic kinematics introduced in the 
medium can enhance the polarization. 

In this Letter, we compute the polarized flux of the emergent Lyct from an optically 
thick and expanding slab. This result is applied to a hemi-spherical shell that is expected 
in various systems including primeval galaxies exhibiting P-Cygni profiles. 

2. Theory 

2.1. Anisotropically Expanding Slab 

In AL98 is given a brief discussion about the line radiative transfer in an expanding 
and optically thick medium (see also Osterbrock 1962, Adams 1972, Rybicki and Hummer 
1978, Sobolev 1960). We describe the model adopted in this work and introduce the 
computation procedures. 

First, we consider a slab expanding in the normal direction, for which the bulk 
velocity law is given by ~V(z) = Hzz where z— direction is chosen to coincide with the 
normal direction and the velocity gradient H is constant. Here, the coordinate system is 
chosen such that at z = the bulk velocity vanishes. We also assume that the medium is 
characterized by a uniform density and a fixed temperature. 
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The scattering optical depth of a given line photon with frequency deviation x in units 
of the Doppler width is given by 



r(s) = J\-f° du e- 2 

COS 9 J-oo 



, { x — u\ i / x — u — s cos 9 
tan" 1 - tan" 1 



a J \ a 

where the parameter s is used for the path length / of the photon by the relation s = Hl/v t h, 
and 9 is the angle between the z— axis and the propagation direction of the photon. Here, 
r = ^fosc^pi^j * s ^ ne Sobolev type scattering optical depth. We fix the velocity gradient 
H, and the bulk velocity difference at the bottom and the top of the slab is given by 



AV = Hd = 100 J kins" 1 , (2) 

where d = s max v t h/H is the slab thickness. Other parameters have their usual meanings 
(see AL98 for detail). 

The polarization state of a given photon is described by a density matrix and the 
circular polarization is zero in this work because of the azimuthal symmetry. The density 
matrix associated with the scattered Lya photon is dependent on the scattering type and in 
particular if the scattering occurs in the damping wing, then the phase matrix is identical 
to that for the Thomson scattering or the classical Rayleigh scattering. A more detailed 
description of the density matrix associated with the scattered photon as a function of the 
incident photon states is provided by Ahn & Lee (1998). 

We locate the incident photon source at the bottom of the slab and perform a Monte 
Carlo simulation to get the emergent flux from the other side of the slab. The initial 
incident profile is chosen to be a Gaussian of width a = 5v t h- The linear degree of 
polarization is computed with respect to the normal direction of the plane spanned by 
the wave vector k of the emergent photon and z-direction. In this work a positive linear 
degree of polarization represents the electric vector lying in that direction and a negative 
polarization perpendicular to it. 
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2.2. Expanding Hemispherical Shell 

The preceding slab analysis is applied to deal with an expanding spherical shell. 
Assuming that the spherical shell can be decomposed into slabs with normal vectors 
distributed spherically, we integrate the emergent flux from each slab component to get the 
profiles and the polarization. The hemispherical geometry is chosen because it describes 
shell structures in many objects including neutral shells surrounding an HII region such as 
NGC 1705 (Meurer et al. 1995), or galactic supershells (e.g. Kamphuis et al. 1991, Heiles 
1979, 1984, Koo 1991, 1992, Reach 1993) possibly obscured by optically thick dust lanes 
(Ichikawa et al. 1994, Scoville et al. 1998). 

3. Results 
3.1. Lya Line Formation and Polarization 

EDITOR: PLACE FIGURE HERE. 

In Fig. 1 we show the main results. We collect those photons emerging with 
/i = cos# = 0.5 ± 0.05. The overall polarization increases as AV increases. The polarized 
flux is significant in the red part when r = 10 4 , s max = 10 and r = 10 5 ,s max = 1. A 
typical degree of polarization is found to be 10 per cent in those cases, and the polarization 
direction is perpendicular to the normal direction of the slab. We note that only a small 
degree of polarization is obtained when tq < 10 3 , where scatterings rarely occur in the 
damping wing regime. 

Line photons with moderate frequency deviation (|x|^4) suffer a large number of 
resonance scatterings at a local site and the radiation field gets isotropized (e.g. Lee 1994). 
After a severe frequency diffusion, scatterings occur dominantly in the wing just before 
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escape. Therefore, diffusion in the real space as well as in frequency space is now important, 
and thereby the escape of the line photons is dependent highly on the detailed geometry 
of the scattering medium. Thus, being characterized by the Rayleigh phase function, the 
wing scattered photons can be significantly polarized. Furthermore, an introduction of an 
anisotropic expansion enhances spatial excursions into the normal direction and leads to an 
increased degree of polarization. 

In a thin electron scattering atmosphere, the emergent photons are polarized in the 
normal direction to the slab plane (see Phillips & Meszaros 1986), and we point out that a 
similar analogy explains the polarization direction of the emergent Lya photons shown in 
Fig. 1, because they suffer a few damping wing scatterings just before escape. 

In the case of semi-infinite electron scattering atmosphere of which a detailed 
investigation was made by Chandrasekhar (1960), the emergent continuum radiation 
is polarized in the slab plane direction, because the electric vector associated with the 
emergent radiation relaxes to the slab plane direction. A similar phenomenon is expected 
in the case where photons suffering a large number of damping wing scatterings just 
before escape are dominantly obtained. However, due to frequency redistribution they are 
distributed in a large frequency range and mixed with other line photons suffering fewer 
wing scatterings. Hence only near the line center a small flux of photons with a large 
number of wing scatterings can be found. This observation leads into a possible polarization 
flip near the line center, which is marginally seen in Fig. 1 in the cases r = 10 4 , s max = 10 
and t = 10 5 , s max = 1. 

A more detailed analysis about the basic physical processes will be discussed in a 
subsequent paper (Ahn and Lee 1998). 
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3.2. Hemispherical Shell 

EDITOR: PLACE FIGURE @ HERE. 

In Fig. 2 we show the averaged degree of polarization expected from a hemispherical 
shell described in section 2.2. Here, the obscuration by an optically thick dust lane is 
assumed to extend from the equatorial region to the region with colatitude 9 = 0°, 30° and 
60°. The adopted paramteres are r = 10 4 , s max = 10. 

Significant polarization and polarized flux are found in the red part and this depends 
on the obscuration fraction. In the case of full obscuration (8 = 0°) we obtain negligibly 
small polarization. The linear degree of polarization is obtained ~ 5 per cent where the 
obscuration fraction is 1/2. Therefore, a positive detection of polarized Lya with P-Cygni 
type profile from a cosmological object can be a strong test of the expanding shell structure 
with obscuration. 

4. Summary and Discussion 

It seems a general consensus that the P-Cygni type Lya emissions are originated from 
expanding envelopes of H II regions, which are indicative of the massive star formation. 
These are often obscured by dust lanes or thick molecular disks (Ichikawa et al. 1994, 
Scoville et al. 1998). 

In this Letter we computed the polarization of the Lya photons that are transferred 
through an optically thick and expanding neutral hydrogen layer. Anisotropic expansion 
and high column density are coupled to enhance scatterings in the damping wing into the 
direction corresponding to the largest velocity gradient, which results in highly polarized 
emergent flux. In particular, in a spherical shell with column density ~ 10 20 cm -2 and 
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expansion velocity ~ 100 km s 1 we find that the averaged degree of polarization of the 
emergent Lya line photons reach as high as 0.05 when /i = 0.5. 

There are three interesting classes of primeval objects showing P-Cygni type profiles; 
the DLA candidate galaxies including DLA 2233+131 (Djorgovski et al. 1996, Lu et al. 
1996, 1997) and DLA 2247-021 (Djorgovski 1997), the Lyman break galaxies at 3 < z < 4 
(Steidel et al. 1996, 1997, Lowenthal et al. 1997), and the remote galaxies observed in the 
gravitational lens surveys (Franx et al. 1997, Frye et al. 1997, Trager et al. 1997). 

Firstly, several DLA galaxies with 3 < z < 4 are listed including DLA 2233+131 with 
z — 3.15 by Djorgovski et al. (1996,1997) who concluded that these objects are progenitors 
of normal disk galaxies today. 

In view of the point that the galactic rotation would erase the P-Cygni structure, the 
Lya emitters or the giant H II regions might be concentrated in a compact region. It is 
noted that they also exhibit the P-Cygni type Lya profiles, which is one of prominent 
characteristics of nearby dwarf starbursts (e.g. Meurer et al. 1995, Kunth et al. 1996, 
Legrand et al. 1997, Leitherer 1997). According to Meurer et al. (1995), there is evidence 
for existence of a dust lane or a thick H I disk component partly obscuring the UV sources 
or the super star clusters. 

Secondly, interesting astronomical objects showing P-Cygni type profiles are the Lyman 
break galaxies at z ~ 3 discovered by Steidel et al. (1996) (see also Lowenthal et al. 1997). 

From ~ 75% of the Lyman break galaxies obtained from their recent survey, Steidel 
et al. (1997) detected Lya emission lines which are weaker than expected from the UV 
continuum luminosities. They ascribed this to the resonant scattering (see Kunth et al. 
1998). 

The third group includes the remote galaxies which can be observed using the 
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gravitational lens effect of the clusters of galaxies (Franx et al. 1997, Frye et al. 1997, 
Trager et al. 1997). Majority of the obtained profiles show the prominent P-Cygni type. 

One other important application may be found in the first generation star clusters 
investigated by Haiman and Loeb (1997a,b). These primeval objects are expected to be 
observed at the redshift of 5 < z < 10 and therefore the Lja emission will be located most 
possibly in the near IR band where the NGST can be used. The neutral hydrogen layers 
enveloping the super star clusters may trace the Hubble-type expansion, and the mechanism 
of the Lya line formation may be studied in a similar way introduced in this work (Ahn & 
Lee 1998). 

We conclude that the Lya lines emerging from an expanding and optically thick 
medium are measurably polarized by up to 5 per cent and therefore a good constraint on 
theoretical models is provided by band polarimetry. 

We thank Dr. Leitherer and Dr. Meurer for kindly providing their preprints and 
reprints. We are very grateful to Prof. Blandford for his invaluable comments about this 
work. 
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Fig. 1. — Linear degree of polarization (points with error bars), emergent flux (thick 
solid lines), and polarized flux (thin solid lines) from linearly expanding slabs with tq = 
10 2 , 10 3 , 10 4 and s max = 1, fO. 



Fig. 2. — Linear degree of polarization (points with error bars), scattered flux (thick solid 
lines), and polarized flux (thin solid lines) from a hemispherical shell that is made of local 
slabs with r = 10 4 and s max = 10. The obscuration by an opaque component extends from 
the equatorial part to the region of colatitude 6 = 0°, 30°, 60°. The viewing angle is chosen 
so that the averaged degree of polarization is maximal. 
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